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One-loop SUSY QCD radiative correction to gb tll~ cross section is calculated in the Minimal 
Supersymmetric Standard Model. We found that SUSY QCD is non-decoupling if the gluino mass 
and the parameter /x, At or A}, are at the same order and get large. The non-decoupling contribution 
can be enhanced by large tan /3 and therefore large corrections to the hadronic production rates at 
the Tevatron and LHC are expected in the large tan (5 limit. The fundamental reason for such 
"^v^ ' non- decoupling behavior is found to be some couplings in the loops being proportional to SUSY 

, mass parameters. 

o 
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I. INTRODUCTION 

Although the Standard Model (SM) is phenomenologically successful, it is arguably an effective theory and new 
physics must exist at high energy scales. Among all elementary particles predicted by the SM, top quark and Higgs 
boson may hold the key to new physics since they are most related to the electroweak symmetry breaking. An 
intensive study of the properties of top quark and Higgs boson will be one of the primary tasks of particle physics in 
^ ^ the new millennium. 

. So far the most intensively studied new physics model is the Minimal Supersymmetric Standard Model (MSSM) [Q. 
This model predicts the existence of five Higgs bosons, iJ", and i/^, all of which couples to top quark. 

Compared to the couplings in the SM, the coupling thH~ is an utterly new coupling. Studies show that this 
coupling is sensitive to quantum corrections and may be a good probe of the MSSM. Although this coupling could 
\ be measured from top quark decay process t — > if the charged Higgs is sufficiently light, the direct production of 
(~| ' a top quark associated with the charged Higgs boson through the subprocess gb tH~ at hadron colliders will be 
a good probe for tbH^ coupling In this work we calculate the one- loop SUSY QCD corrections to this process 
^ owing to the following motivations. Firstly, if the charged Higgs boson is heavy, 171^+ > 'nit + as a main charged 
^ Higgs production channel the process gb tH~ will provide sizable cross section at the Tevatron and LHC. 
f~| : The supersymmetric radiative corrections, especially the SUSY-QCD corrections, to this high energy process may 
be significant, as was found for other similar processes §-|l§. Secondly, some recent studies [|8|-|l0[ showed that 
• ^ , SUSY-QCD may be non-decoupling in some processes involving Higgs bosons. As is well known, the decoupling 
■ theorem ||ll| states that under certain conditions in a given quantum field theory with light and heavy particles, if 
5h \ the heavy particles are integrated out to all orders in perturbation theory, the remaining effective action to be valid 
. . . 1 at energies much lower than the heavy particle masses does not show any trace of these heavy particles. If SUSY 
QCD is non-decoupling in some cases, we need a proper understanding and thus we need to further investigate such 
non-decoupling property of SUSY QCD. gb — > tH~ will be an ideal process for this purpose. 

This paper is organized as follows. In Section || we present the formula for the one-loop SUSY-QCD corrections 



> 



to the gh — + tH^ process. In Section IH we scan the parameter space of MSSM to estimate the size of SUSY-QCD 
corrections. In Section IV we study the decoupling behavior of SUSY QCD. A discussion on how the decoupling and 
non-decoupling take place is also given. Finally, the conclusions are summarized in Section ^. 

II. CALCULATIONS 

The subprocess gb — > tH^ occurs through both s-channel and t-channel. The tree-level amplitude is given by 

Mo^M^^'Hm^'\ (2.1) 

where M^^^ and M^p represent the amplitudes arising from the s-channel diagram shown in Fig. 0(a) and the t-channel 
diagram shown in Fig. 0(6), respectively. Their amplitudes can be expressed as 
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where Pr^l = (1 i 75)/2, and pt, Ph and k are the momenta of the outgomg top quark, the incoming bottom quark 
and the incoming gluon, respectively, s and t are the subprocess Mandelstam variables defined by s = {pb + fc)^ = 
{Pt + Ph-Y ^'^d i = [pt — ky = {ph- — PbY- T"- are the SU{?>) color matrices and tan/3 — V2/vx is the ratio 
of the vacuum expectation values of the two Higgs doublets. The constants rj^^t are defined by rji, — rrib tan [3 and 
r]t = mt cot (3. 

The one-loop Feynman diagrams of SUSY QCD corrections are shown in Fig. |l|(c)-(o). In our calculations we use 
dimensional regularization to control all the ultraviolet divergences in the virtual loop corrections and we adopt the 
on-mass-shell renormalization scheme. The renormalization condition for the coupling constant Qs is similar to that 
for the coupling constant e in QED, i.e., the coupling of the photon (gluon) to a pair of fermions is required to recover 
the tree-level result in the limit of zero momentum transfer. This condition yields 
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(2.4) 

Sgs and = ^/^A^ with 5" denoting 



where 5gs and Z| are the renormalization constants defined hy g1 = ^ 
bare coupling constant and A^^ bare gluon fields (color index suppressed). 

Including the one-loop SUSY QCD corrections, the renormalized amplitude for gb — > tH~^ can be written as 
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where SM represents the one-loop SUSY QCD corrections given by 



(2.5) 
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Here (5M^i(*\ (5M^2(s) g^j^j SM'^^^^ represent the renormalized vertex gbb, tbH^ and the renormalized propagator in 
the s-channel diagram, respectively. Similar definitions are for (5M^'^(*\ JM^^^*) and SA'P^^^ in the t-channel diagram. 
The contribution of the box diagram is denoted by 5M^°^ . Each 5M^ can be decomposed as 
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(2.7) 



where the coefficients C' and the form factors F^^ are given explicitly in Appendix A and B, respectively . We have 
checked that all the ultraviolet divergences canceled as a result of renormalizability of MSSM. 
The amplitude squared is given by 
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Here the color factor Nc — 3, and hh and hn can be found in Appendix A. 
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FIG. 1. Foynman diagrams of gh tH^ with one-loop SUSY QCD corrections: (a) and (6) arc tree level diagrams ; (c) — (e) 
are one-loop vertex diagrams for s-channel ; (/) — {h) are one-loop vertex diagrams for t-channel ; (i) — (m) are self-energy 
diagrams; (n, o) are the box diagrams. 



The cross section for the parton process gb tH is 



(2.11) 



with 
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i + - J ± ^[s - [mt + tuh- Y] [s - {mt - mn- )^] | • (2-12) 

The total hadronic cross section for ppipi pp) — > tH^ + X can be obtained by folding the subprocess cross section 
a with the parton luminosity 

cr(s) = / dT — a{s^sT), (2.13) 
Jro dr 

where tq = (mt + rriH-Y /s, and s is the ppipr pp) center-of-mass energy squared. dL/dr is the parton luminosity 
given by 

dT dr 

^ = I -[f^{x,Q)f!{r/x,Q) + (5 ^ fe)], (2.14) 

where and fP are the bottom quark and gluon distribution functions in a proton, respectively. In our numerical 
calculation, we use the CTEQ5L parton distribution functions |]l2| with Q = nit + mu- ■ 
To show the size of the corrections, we define the relative quantity as 



where ctq is the tree-level cross section. 



^SQCD = - — (2-15) 
Co 



III. NUMERICAL RESULTS 



Before performing numerical calculations, we take a look at the relevant parameters involved. 

For the SM parameters, we took niw = 80.448 GeV, mz = 91.187 GeV, mt = 176 GeV, mt = 4.5 GeV, and used 
the two- loop running coupling constant as{Q). 

For the SUSY parameters, apart from the charged Higgs mass, gluino mass and tan/3, the mass parameters of stops 
and sbottoms are involved. The mass square matrices of stop and sbottoms take the form {q = t or b) u3 
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where 



m?^ = mQ + rn^ — ^2(2 + ^^"^ ^w) cos(2/3), 

TO?^ = m^ Q + -I- eqm\ sin^ Ow cos(2/9), 

X ^ f Af - ^ cot /3, for q^t, ,g 2) 

— ^tan/3, for g = 6. ^' ' 

Here m^, mjj and m^^ are soft-breaking mass terms for left-handed squark doublet Q, right-handed up squark U and 

down squark D, respectively. Ab{At) is the coefficient of the trilinear term HiQD {H2QU ) in soft-breaking terms 
and /i the bilinear coupling of the two Higgs doublet in the superpotential. Thus the SUSY parameters involved in 
stop and sbottom mass matrices are 

mg , , TO^ , , Afc , ^, tan /3. 

The mass square matrices are diagonalized by unitary rotations {q = t or b) 

cos 9g sin Of 
9q cos 9q 



W= ( ^^r-^^ '^'"'"^ ] (3.3) 
> — sm On cos d ' ^ ' 



which relates the weak-eigenstates {qL,qR) to the mass eigenstates (gi,g2)- Then the stop and sbottom masses as 
well as the mixing angles are obtained by 
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To determinate the mixing angles completely, we adopt the convention in |14[ which sets 9q = 7r/4 if rriq 
shifts 7r/2 to 9q if ruq^ > rriq^. Thus 6q lies in the range —j < 9q < jTt. 

To find out the size of the one-loop SUSY QCD effects, we performed a scan over the nine-dimensional parameter 
space: TOq, rrifj, rrif), At, Ab, /U, tan/S, rriH-, rUg. In scanning we restricted mu-, At, A^ and ^ to the sub-TeV 
region and required tuh- > 150 GeV. Other mass parameters are assumed to be smaller than 5 TeV. In addition, we 
consider the following experimental constraints: 

(1) > and a large tan/3 in the range 5 < tan/3 < 50, which might be favored by the recent muon g 
measurement | |T5| ]. 

(2) The LEP and CDF lower mass bounds on gluino, stop and sbottom [|l^ 

> 86.4 GeV, > 75.0 GeV, m-g > 190 GeV. 

The scan results are plotted in the plane of Asqcd versus 9b in Fig. |[ 
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FIG. 2. The scatter plot in the plane of Asqcd versus 9b- The scan was performed over nine SUSY parameters. 9b is in 
unit of TT. 



From Fig. |^ one can see that in most parameter space the mixing of sbottoms is small while the one-loop SUSY 
QCD effect can be quite large. In some part of the parameter space, the correction size can be larger than 20% which 
cannot be neglected in the study of this process at LHC. 
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IV. DECOUPLING PROPERTY OF SUSY QCD 



To find out if SUSY QCD is decoupling from the process gb tH^ in tlie large limit of SUSY mass parameters, 
we fix the charged Higgs mass as rriH- = 250 GeV and consider the following scenarios. 

(1) Scenario A: All squark (collectively denoted by ms), gluino masses and fj, or A parameters arc of the same 
size and much heavier than the electroweak scale, i.e.. 



ms 



fi or At or A^ ^ ms » Mew 



(4.1) 



In this case, both mixings in the sbottom and stop sectors reach their maximal values, i.e., 9t ~ ±-|, Ob ~ 
As shown in Eqs. ( |6.3 -6.5), the couplings aij in the vertex H~iibj are proportional to the linear combination of 
fi + Ab tan l3 and fi + At cot (3. Considering that the couplings are proportional t o ms as /i or At^t gets large as 
ms, and the loop scalar integral functions Co goes to — 1/2to| as ms ^ s (see Eq. ( 3.18 ) and Eq. (B8) in Ref. [^), 
one can infer that the terms Uij Af^^'^^ CQmg which arise from the vertex correction to H~th do not vanish but go to 
a non-zero constant, showing a clear indication of non- decoupling behavior. In fact, such non-decoupling behavior will 
happen as long as the gluino mass and or A parameters are of the same order, not necessarily degenerate. Actually, 
from the expansions of the three- and four-point loop integrals in the asymptotic large mass limit in Appendix C, one 
can see this fact. 
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FIG. 3. Non-decoupling behavior of Asqcd with uiq = rrijj = m^, = rrig = Ab = At — jJ, = ms and for different values 
of tan/3. Corrections at the Tevatron with ^/s = 2 TeV (solid lines) and at the LHC x/s = 14 TeV (dashed lines) are plotted 
respectively. 



As illustrative examples, we plot the dependence of the SUSY-QCD correction to gb — > tH~ on the common SUSY 
parameter ms in Fig. |^. From this figure one can see that the non-decoupling behavior indeed happens. As for the 
dependence of the non-decoupling effects on tan/3, it is quite involved and complicated. For the parameter values 
chosen in our numerical examples, the corrections are enhanced by tan/3. 

(2) Scenario B: The gluino mass and /x are of the same order (collectively denoted by ms) and get much larger 
than squark masses and electroweak scale, i.e.. 
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rrig ^ ^:$> ms 
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To keep stop and sbottom masses from getting large, we can set At ~ /icot/3, Af, ~ /itan/3. In this scenario, no 
mixings occur in the sbottom and stop sectors. Apart from the vertex correction discussed in Scenario A, the terms 
such as aijA^j^'^^ DgTrig in the box contribution (see Eq. (3.12) and (3.24)) do not vanish either since the four-point 
integral functions D^^ 1/"^^ {I ~ Ij 2, 3) when m? ^ s. So in this case the SUSY QCD is non-decoupling. This 

differs from Scenario A where -0(0,1;) 

Although the non-decoupling effects can arise from more diagrams in this case, the reason is the same as in Scenario 
A, i.e., the couplings H~tibj are proportional to /x. To prove this point fully, now let's focus our attention on the terms 
arising from the corrections to the ghh and ggg vertices which are also likely to give the contributions to non-decoupling 
effects in both scenarios we discussed. Firstly, for example we consider the term 



377* [mlCl - 2CI4 + 1/2] 



(4.3) 



m fI'^'\ Fr om Eq. (|6.13|), (|6.16|),(|6.18|) (|6.19|) (|6.2l|) and (|6.22|), we draw a conclusion that the term indeed cancels 



out. Further, it is easy to find that all terms related with Af^^ are zero in the asymptotic large mass limit. This is 
also valid in the Scenario C and Scenario D we will study. Secondly, we exam the terms such as aijA^j^'^^ Dgm^ in 
Eq. (6.12). The asymptotic form of the function Dq in large mass limit is always proportional to -\ and different 

from _Dq which is proportional to in the Scenario B. So the terms don't cause non-decoupling either. 
(3) Scenario C: Only the gluino mass is very large than other SUSY parameters and the electroweak scale 
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FIG. 4. Behavior of Asqcd in the large mg limit with fixed mg = rufj = mf, = Ai, = At = fJ, — 1 TeV and for different 
values of tan (3. The solid and dashed lines correspond to corrections at the Tevatron with y/s = 2TeV and at the LHC 
^/s = UTeV , respectively. 



In this scenario, to simplify the calculation we assumed ms = rUg = m^y = = fi = At = Af, — 1 TeV. As 
shown in Fig. g the SUSY-QCD decouples. The reason is the couplings in the vertex H Ubj arc fixed, and in 
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this case the scalar function Hke Cq is proportional to ^ log ^ 1 + when TO? > s (see Eq. (|6l^ , |6!26|) ), thus 

9 g L g J 

2 

the SUSY QCD correction goes hke — log 

^ 3 

In the scenario, since s can be up to the large collider beam energy squared s, therefore, not only the logarithmic 
dependence on the large mass parameter but also the large collider beam energies are responsible for the slow decou- 
pling of the gluino, especially, at LHC, as shown in the figure. This is different from the previous studies in various 

decays 

(4) Scenario D: Only squark masses are of the same order (collectively denoted by ms) and very large than other 
SUSY parameters and the electroweak scale 

As shown in Fig. ||, the SUSY QCD also decouples. In this case, it decouples much faster than in Scenario C where 
only gluino mass gets large. This can be understood easily because in this case the couplings aij in the vertex H~tibj 



and the mass of gluino are both fixed, so that the SUSY QCD correction goes like ~ Cli 
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FIG. 5. Behavior of Asqcd in the large squark masses limit with fixed mg = Ai, — At — fi = 1 TeV and for difi'erent values 
of tan /3. The solid and dashed lines correspond to corrections at the Tevatron with ^/s = 2TeV and at the LHC ^/s = 14TeV , 
respectively. 



Some remarks on the SUSY QCD correction to the process of gb — > tH^ are in order: 

1. From the above analysis we find that the fundamental reason for such non-decoupling behavior of SUSY QCD 
in the process gb tH^ is that some couplings like H~iibj are proportional to SUSY mass parameters. This is 
similar to the non-decoupling property of the heavy top quark in the SM, where the top quark Yukawa couplings 
are proportional to top quark mass. 

2. The non-decoupling behavior shown in Fig. |^ is in agreement with previous studies of SUSY-QCD corrections in 
some decay process |^,^. In particular, the correction Asqcd shown in this figure as a function of the common 
scale ms and tan/3 looks quite similar to the corresponding corrections in the partial decay width T{H~^ — s- tb), 
as given in Fig. 2 of Ref. The same kind of similarities are also found between Fig.|| and Fig. 6 of Ref. [||, 
and between Fig.^ and Fig. 7 of Ref. Although the processes are different, the fundamental reason for such 
non-decoupling behavior is the same. 
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3. From Figs. |^|| one sees that the size of SUSY-QCD corrections can be quite large for large tan/3. Note that 
when one-loop effects are too large, higher-level loops must be also calculated. We refer the reader to Ref. jl^] 
where some techniques of resummation for a better convergence are proposed. 

4. As shown in Ref. jH, the genuine QCD corrections to this process is also sizable, which can enhance the 
production rate by 40% ~ 80% when charged Higgs mass and tan/3 vary in the range 180 ^ 1000 GeV and 
2 ^ 50, respectively. It is clear that the SUSY-QCD correction evaluated in this work is comparable in size to 
the genuine QCD corrections. It is noticeable that the genuine QCD corrections are always positive whereas 
the SUSY-QCD corrections are negative in most SUSY parameter space. 



V. CONCLUSIONS 



In this work we have evaluated SUSY QCD radiative corrections to the gb tH^ at Tcvatron and LHC. We have 
found that in some parameter space the one-loop SUSY-QCD correction can be quite large and cannot be neglected. 
We have discussed in detail on the decoupling behavior of the corrections in the large SUSY mass limit, and found that 
with fixed gluino mass the one-loop SUSY QCD corrections decouple; while non-decoupling occurs when gluino mass 
and ^ ox A parameters both get large. The non-decoupling behavior of the SUSY-QCD corrections in the process 
gb — > tH~ is similar to the ones found in the literatures for the Higgs particles and top quark decays. We pointed out 
that such non-decoupling behavior arises from the H~tibj vertices which are proportional to SUSY mass parameters, 
as stated in some previous literatures [p|-p^. Such large non-decoupling effects may play an important role in the 
indirect search for SUSY from the production of a top quark associated with a charged Higgs boson at Tevatron and 
LHC. 
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VI. APPENDIX A 



Here we hst the coefficients C', scalar functions /in' and the vertex V{H iibj) — igaij / V^niw — 1,2) needed 
in our calculations. 



Coefficients C 



1 1 



, C^°=^ = 1. 



(6.1) 



Scalar functions h 



(0 



Amtr]t{2pb ■ k -p(') ■ pb) - 4mbi]b{p'^^^ ■ pt +Pt ■ k), 

2T]t{2pb ■ kpb ■ Pt - mlpt ■ k - 2p^^^ ■ pbPb ■ Pt) + 2mbmtr]b{pb ■ k - 2p('^ • pb), 
h\}' = 2r]t{m^pb ■ k - 2p('' • ptPb ■ Pt) + 2mbmt7]b{pt ■ k ~ 2p'^^^ ■ pt), 

4?7t(p*^'^ • PbPt ■ k - p'^'^ • kpb -pt-Pb- kp'^'-'' ■ Pt - 2pb ■ kpt ■ k) - AmbmtrjbP^^^ ■ k, 
h'^> = Amt-qtPb ■ k{pb ■ k - p'-'-^ ■ pb) - 'inibrjbp^^'^ ■ PbPt ■ k, 
"^mtritPb ■ k{pt ■ k - p'-^^ ■ pt) - AnibribPt ■ kp'^^^ ■ pt, 

"•2,4,6,8,10,12 — "-1,3,5,7,9,11^/6 

where the index I represents the two channels s and t, and p*-*-* = pb, p^*^ = pt- 



(6.2) 
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• Couplings of H iibj 
The H~iibj interaction terms can be parametrized s 

^H-kh = ^^^ Q^'J i^^j = 1' 2), (6.3) 

where 

(^ij = ^iiI^ji9LL + R*aR%9RR + R*aR%9LR + R*^R]i9rl (6.4) 

with 

gLL = —m^ sin 2/3 + ml tan/3 + cot /3, 
ffflij = msmt (tan /3 + cot/3), 
gLR = mb{iJ. + Ab tan /3) , 

ffijL =mt{ii + At cot /3). (6.5) 



APPENDIX B 

The form factors raise from the renormalized vertices and propagators of s-channel and t-channel, as well as 

box diagram given as following. 

• The renormalized vortices of s-channel: 

Fi^^(^) = -3776 {Almgpk ■ kC] + Alm^ [sC^ + ml{Cl + 4Ci\ + 2Cl,) 

+2pb ■ fc(Ci\ + + 2CI3) - 2m?Ci - 1 + + rribPb ■ k{Cl + Cl^)] } 



2 



+^^^ir?6m6(7|4 - ^T]tmb{5Z^ji - 5Zl), 
<J 9s 

fI"^'^ = \nt {- K(Ci\ + Cl^ • k(C\^ + Ci3) + + Ai,m^m5(Ci + 

+2^2^ {Cl^+Vb ■ k{Cl, + C2I3)] } - 3ryt { [m2(C2 + 2C2i + Cl) + m^Cl - 2Cl^ + 1/2] 

+2Almbmg{Cl + C^^) + 2Al [ml{Cl + 2C?i + C|i) - mfC^ - 1/2 + 2Cj,] ] 



2 



+ ^^m25Z^R, 

9s 

F^'^'^ = \m{Cl^ - '2AlCk) - Ivt {mliCl - C^) + 2pb ■ k{C!, + Cl,) - m\Cl + 20%^ - 1/2 

-2Almbm-,Cl - 2Al [ml{C', + 2Cf, + C|i) + 2pb ■ k{C!^ + Cl,) - m\Cl - 1/2 + 2C|J } 

16^2 



-2-^*^^fl' 
i/s 



^^J"'^'^ = \nh [-m6(Ci\ + Cli) + A'umg{Cl + C^^) + 2A2,m6(2Cii2 + 2C\^ - - C\^)\ 
-3% {mbiCl^ + C|i) + Aj,m§Ci^ + 2A2.m6(Ci^ + C|i - 2Ci22 - 2CI3)] , 

pV.is) ^ pMs) ^ ^2^ ^ _^2^^ ^ (g_g^ 

where the Feyman integrals are defined as = C {pb,k,mi ,mg,mg) , = C [—pb,—k,mg,mi^,,mi^,), = 
C [-Pt, -Ph- , m-g, mi. , mg^ ) and 
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Here Al^ ^,^^ = (-l)'sin26lt,6, Af^ ^^^. = cos2(9t,6. Ay are defined as A}[ = 2i?|2'Rj2, ^if = '^RliRji, 

-2Rl2Rj, and = -2RlR%.' 

• The renormalized vertices of t-cliannel: 



r^i(t) 



-AlrjtmtC^^ + irit {mtPt ■ k{C^ + Cfj + Al^mgPt ■ kCl 

-Almt [iCi + m\{Cl + 4Ci\ + 20%^ - 2pt ■ k{C!, + 2Cl^ + 2Cl^) - 2mlCl - 1 + Wl^)] ] 

-l^W {6Z'n-5Zl), 
9s 

\m {- K(Ci\ + Ct^) - Pt ■ k{Ct2 + Ct^) + Cl^] + A\,mgmt{Ct + Cfj 

+2AI [-Cl, +pt ■ k{Ct2 + C^s)] } - 3% {m^Cl + 2Cl, + C|i) + mpl - 20%^ + 1/2 

^2A\^mtmg{Gl + Cfj - 2Al {m\{Cl + 2C\^ + - m?Co^ - 1/2 + 20%^ } 

\^t{G\^ + 2A2.C4J + ^^2(c^5 _ (.5^) + • k{C\2 + CI3) + mlGl - 20%^ + 1/2 

+2A\^mtmgCl - 2Al {m\{Cl + 2C\^ + - 2^* • kifll^ + CI3) - m\Cl - 1/2 + 2CI4] } 
167r2 



+ 



pVi(t) _ 

2,6,8 ~ 
^11 - 



„2 ^f^^l?' 
as 

V'i(t) 



i% [mt(Ci\ + C^i) - Alm^,{Cl + Cfi) - 2Almt{2C\2 + 2CI3 - Cfi - C\^\ 
+3r;t [mt(Ci\ + C\^) + A,\m§Ci\ + A'i2mt(Ci\ + C|i - 2C7f2 - 2CI3)] , 



-'^12 — -^11 \J)t Vb, —>■ -^bi) , 

wliere the Feyman integrals are defined as 
C [-Pb-,PH-,mg,mi,,mi^ and 



C {-pt,k,mg,mt.,mi.), = C {-pt,k,mt.,mg,mg), 



?V2(t) 



F, 



V2(t) 



7V2(t) 



Fr 



V2{t) 



——aijAjj^Pt ■ kC 12 
F^'^''> {R^L), 

|ay [A][mb{Cf, - Cf^) + Al^mtCf^ - A^fm-gCt] 



167r2 



-Uij [Alj-mbiCt, - Cf^) + A\fmtCl2 - A^^m-gCl] + ^ 



+ 5ZI + 2 



6mt 
mt 



:SZi 



:SZ'r 



5mt 
mt 



F-f=F-f {R^L,r,t 



Vb, 



5m,t Smb 



mt 



mb 



11 



The renormalized propagator of s-channel: 

167r2 



■,is) 



Fr 



.s{s} 



2,4,1 



9 s 
167r2 

—2-'nt 



9 s 

16^2 



5: 



-^1,3,7 ^ rib,R^ L) . 



(6.10) 



• The renormaUzed propagator of t-channel: 

167r2 



(t) 



5f 

16^2 



-'7t 



9t 
167r2 



2TOtpt • k {±1 + 
2fe + m2s*^) +Amm 



2 



(6.11) 



• Box diagram contribution: 



jpbox 



rpbox 
^5 



TT^box 
^7 



TT^box 
^9 



box 



T^OC 
TTihoX 

-^2,4,6,8,10,12 



1 



IR 



^27 + 2 (™?^22 + - Pt • kDl^ +{pfk^pb- k)Dl^ 



+ {pfPb+Pfk- m^t)Dle - pt ■ k{D\^ - D\^) - -mlol 



-D. 



27 



ij 

1 



2-D25) 



-a 



AlfmtiDt, + Di,~Di,) 
{9 [Ai/m,(i5i3 - Dl,) - A\^mt{Dl^ + Dl^ - 2Dl,) + Ajj^m.g{D\^ 



+A^fm^,{Dl + Dl,-Dl^)], 



+A\fm,{Dl^ - Dl,) ~ Afm,Dl^ + Alfm^Dl^] 



[A\j^m,D\^ + Alfi7i,{D\^ - D\^) - Ajfm^.Dl] , 



•'24 



-a,,A\^' {-%D\^ - Dl,) + Dl, - Dl, - D 

pbox 



Die} , 



-f'l,°3!5,7,9,ll ^ F) 



(6.12) 



as 



where four-point functions are defined 

= D (je, —pt, —pH- ,n^g,n^g,'m'iiT''^b^ ^ ^ D (^—pb, —k,pij- ,mg,mf^ ,mf^ ,mj:.^ . Note in the above formula 
we take the convention that repeated indices are summed over. 

All other form factors F^ not listed above vanish. Note that the contributions from diagrams Fig. |^(i),0(j) and |l|(k) 
just give the renormalization constant 6Z2, which is canceled out by 6gs due to the renormalization condition in Eq. 

& . . 

The renormalization constants appearing in the above form factors are given by 



16772 



2mt 



dBi 
dpi 



2 2 2 
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5Zl 



5Z% 



6mtlmt 



167r2 
16 



nib 



2mb 



dBi 
dpi 



rrigA 



1 5i?o 



Pl=rnl + (1 - 2A2JBi I [pI 



2 2 



:^(i?i-^Ai,i?o)(p?,m?,m?V 



167r2 



2 2 2 



(6.13) 



where the color factor Cp — 4/3. The renormahzed self energy contributions from quarks as follows 



yt 



yb 



Idn-^ rrit * 
glCp m 



yt _ 



167r2 
167r2 



D I 2 2 2 

Bq I m^,mg,m~^^ 



il-2Al)BJi,mlml 



-Aid, 
'AhA. 



Bo {s,ml,ml^ 



SZ 



{l-2Al,)B, (s,m?,m|) 



(6.14) 



APPENDIX C 



In this appendix we give the expansions of the scalar loop integrals in the as ymp totic large mass limit. The 
definitions and conventions of the Feynman loop integral functions can be found in |18| . The integrals are performed 
in 4 — e dimension and the divergent contributions are regularized by A = (2/e) — Je + log(47r) — log{ml/ ^q) with 
m/i being the corresponding mass of the heavy particle in the loops and ^q, the regularization scale. 

Under the assumption of ml — max{mD 3> Pi ■ pj = 1,2; k = 1,2,3), we consider special cases used in our 
calculations. For the two-point function Bi{p'^ ,m\,ml), we obtain 



Bi 



A 



25- 



3-<5 



2(1 -J) (l-<5) 



log4+0^^ 



(6.15) 



where 5 — ml /ml. Therefore, the asymptotic form can be expressed as 



B, 



2 ^ 


^ 12m 


i : 


6- 


^1, 




4 + 


" 2 
P 


S- 


-0, 


2 ^ 










4 + 


P 


6- 


-> oo 


2 ^ 









(6.16) 



Note in the case S — s- oo, m^ in Eq. ( |6.16| ) should be replaced by mi. 

For the three-point fmiction C(o,24) (PI7P2, "^i, "^21 "^3) expand them as follows: 



• Case A: m^ — mi, ml = m^ = m§ ± A^ 



Co = 



C2 



1 



1 



1-6 (l-(5)2 
2.-5 



In (5 



25+1 



4(1 -,5)3 2(1 -,5)' 



+ — ^ In 5 - 



[1-5) 4(1-5)2 



24m| 



2 5<5 - 52 



+ 



In (5 



6<5 



(1-5)3 



■In 5 



A2 
+ 



+ 
A2 

ru 
9 



(6.17) 



where p = pi + J32. In asymptotic large mass limit we obtain 



13 



Co 



1 



V ™h / 



12m^ 



1, 

0. 



(6.18) 




■ 1, 

■ 0. 



Case B: m; = mi, m/i = m2 = ma 



Co — 



1 



Co4 = A 



1-S 
-35 



8{1-S) 4(1 -,5)2 
The asymptotic form is given by 



\n5 + 



1 - 5,5 - 2,52 



6^2 



{l-5f 



{l-5y 



+ 



4 

+ 



Co 



1 

"2^ 
1 



-0 



1 + /^ 



6 
5 



(6.19) 



(6.20) 



(6.21) 



C24 = 



mt 



For the four-point functions -Do,ife(pi)P2jP3j m2, to|, m§, TO4) with the assumption of .^^ 
max{pi ■ pj) (i, j = 1, 2, 3; A; = 1, 2, 3, 4), we consider two special cases used in our calculations 

• Case A: ruh = mi, mi = m^ 



(6.22) 



max{m^) ^ = 





1 


Do = 




mlmf 




1 


Dn = 






m2m2 



1 + ^ 



■ + 



■In^ + O 



(4)+o(4 



. [2 + 5,5-^2 35 , , 

>2{^(T37F^(T^^^' + '' 



A2, 



mt 



+ 



mi 



where — max{\m2 — TO3I, |m| — m4|, |to| — mil). In asymptotic large mass limit we obtain 



^ i(6'8'12'24)7^-'-+^(^) 
>>2' 2' 3' 6^ "»l™f I^J- i- >^ V™!/ 



•1, 

■ 0. 



(6.23) 



(6.24) 



• Case B: mh = m\ = m2, m2 = m\ 
Do 

Di3 



1 



{1-sy 



9 + 36 ^2 
+ 



1 

^ \ 4(1 -5)3 
1 f 5 + 5 



(1-5)3 
25-2 



, A2 

ln5 + ( ^ 



2(1 - 5)4 
1 + 25 



A2 



mi J VK 



P 



+ ■ 



2(1-5)3 (1-5)4 



ln5 + 



+ 



(6.25) 
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In asymptotic large mass limit we obtain 



D. 



(0,11,12,13) 




(6.26) 
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